Excitonic condensation of massless fermions in graphene bilayers 
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Graphene, a single sheet of graphite with honeycomb lattice structure, has massless carriers with 
tunable density and polarity. We investigate the ground state phase diagram of two graphene sheets 
(embedded in a dielectric) separated by distance d where the top layer has electrons and the bottom 
layer has holes, using mean-field theory. We find that a uniform excitonic condensate occurs over a 
large range of carrier densities and is weakly dependent on the relative orientation of the two sheets. 
We obtain the excitonic gap, quasiparticle energy and the density of states. We show that both, the 
condensate phase stiffness and the mass of the excitons, with massless particles as constituents, vary 
as the square-root of the carrier density, and predict that the condensate will not undergo Wigner 
crystallization. 

PACS numbers: 



Introduction: Over the past three years, graphene has 
emerged as the unique candidate that provides a real- 
ization of two-dimensional massless fermions whose car- 
rier density and polarity are tunable by an external gate 
voltage Subsequent experimental and theoretical in- 
vestigations have led to a thorough re-examination of 
some of the properties of linearly dispersing massless 
fermions 0, Hf. The truly two dimensional (2D) nature 
of graphene permits control and observation of local car- 
rier density and properties @|. In graphene bilayers, 
the ability to change the carrier polarity of an individ- 
ual layer implies that the interlayer Coulomb interaction 
can be tuned from repulsive to attractive. This raises 
the possibility of formation of electron-hole bound states 
or indirect excitons, albeit with massless fermions as its 
constituents. Properties of such bound states of massless 
particles are an open question; the only other example, 
to our knowledge, is the proposed color superconductiv- 
ity in dense quark matter ■ Graphene bilayers provide 
an ideal and unique candidate for straightforward exper- 
imental investigations of such phenomena. 

A uniform Bose-Einstein condensate of excitons in 
electron-hole bilayers occurs when the interlayer dis- 
tance is comparable to the distance between the parti- 
cles within each layer [§], Q . These excitons have mass 
m ex = m e + to/i where m e (wi/,) is the band mass 
of the electron (hole). At high densities, dipolar re- 
pulsion between the excitons can lead to a condensate 
ground state with broken translational symmetry: a su- 
persolid [lOj. Biased bilayer quantum Hall systems near 
total filling factor one have shown uniform excitonic con- 
densation In this case, the exciton mass is deter- 
mined solely by interlayer Coulomb interaction and is in- 
dependent of the bias voltage 12, 13 1 . These observations 
raise the questions: What is the mass of an exciton with 
massless constituents? Will such an excitonic condensate 
lead to a supersolid if the dipolar repulsion between such 
excitons (with a nonzero mass) is increased? 

In this paper, we investigate the excitonic condensa- 



tion in two graphene sheets embedded in a dielectric and 
separated by a distance d » a (a=1.4A is the honey- 
comb lattice size) so that the tunneling between the lay- 
ers is negligible, but interlayer Coulomb interaction is 
not. The layers have opposite polarity and equal density 
of carriers U2D- We remind the Reader that in graphene, 
in the continuum limit, the length-scale 1/kp and the 
energy-scale Ep are both set by the density of carriers 
1T-2D (kp — ^/im2D is the Fermi momentum, Ep = hvakF 
is the Fermi energy, and vq ~ c/300 is the speed of 
massless carriers). Therefore, the ground state phase 
diagram depends only on one dimensionless parameter 
kpd. This is markedly different from conventional bi- 
layer systems parameterized by (d/as,r s ) where as is 
the band Bohr radius and r s = 1/ ' \f^a? B n2D @], as well 
as biased bilayer quantum Hall systems, parameterized 
by (oI/Ib, A^) where Ib is the magnetic length and Av is 
the filling factor imbalance [HI, 0, EH, [l|| . 

We use the mean-field theory to obtain the ground- 
state phase diagram as a function of kpd. We find that 
a) excitonic condensation occurs at all densities as long 
as kpd ~ 1. b) the condensate properties are weakly sen- 
sitive to the relative orientation of the two sheets (stack- 
ing), c) the superfluid phase stiffness p s and the exciton 
mass have a ^Jtv^d dependence, d) the excitonic conden- 
sate does not undergo Wigner crystallization in spite of 
dipolar repulsion between excitons with a nonzero mass. 

The plan of the paper is as follows. In the next section, 
we present the mean-field Hamiltonian 17[ and briefly 
sketch the outline of our calculations. In the subsequent 
section, we show the results for the excitonic gap A^, the 
quasiparticle energy E^, and the quasiparticle density of 
states D(E). We discuss the density dependence of the 
superfluid stiffness p s and the mass of the excitons. In 
the last section, we show that these results are equivalent 
to absence of Wigner crystallization, and mention the 
implications of our results to experiments. 

Mean-field Model: We consider two graphene sheets em- 
bedded in a dielectric separated by distance d with chem- 
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ical potentials in the two layers adjusted so that the top 
layer (denoted by pseudospin t — +1) has electrons and 
the bottom layer (denoted by pseudospin r = — 1) has 
holes with the same density. We consider two stackings: 
the Bernal stacking that occurs naturally in graphite, and 
the hexagonal stacking in which each sublattice (A and 
B) in one layer is on top of the corresponding sublattice 
in the other layer. Since the Hamiltonian in the con- 
tinuum description is SU(4) symmetric in the spin and 
valley indices, we ignore those indices for simplicity. In 
the continuum limit, the single-particle Hamiltonian for 
carriers in layer r is [l8[ 



H = Skc t (a?iWGfc)c| CQT CkQT 



(1) 



where k is the momentum measured from the K-point 
and a = ± denote the conduction and valance bands 
that result from diagonalizing the Hamiltonian in the 



sublattice-basis. 



(ckctr) is creation (annihilation) 



operator for an electron in band a in layer r with mo- 
mentum k. We point out that for the hexagonal stack- 
ing; c Lr = i c kA T + ae " l9kc Lrl/^ is independent of 
the layer index t. For the Bernal stacking, the creation 
operators in the two layers are related by complex con- 



T ]/V% where 6>k 



jugation, c kar = [cj cAr + ae 
tan (k y /k x ). The interaction Hamiltonian consists of 
intralayer Coulomb repulsion Va(q) = 27re 2 /eQ r and in- 
terlayer Coulomb attraction Ve(q) = — exp(— qd). 

(e is the dielectric constant). Using standard mean- 
field techniques 17|, we obtain the following mean- field 
Hamiltonian 



H = 



k 



£k - A* 
A* 



A k 

£k + 



ek 

;,t 



(2) 



Here el — 



creates an electron in the conduction 



and h . = 



Ck- 



band (a = +) in the top layer (r = +1) 
creates a hole in the valance band (a = — ) in the bot- 
tom layer (t = — ). The term ek contains single-particle 
energy, capacitive Hartree self-energy and the intralayer 
exchange self-energy. The off-diagonal term Ak is pro- 
portional to the excitonic condensate order parameter 
(/i_kek). The eigenvalues of the mean-field Hamiltonian 
are given by ±E k = ±Vte^ - fi) 2 
mean-field states with a real Ak : 



We consider 



A k , and spatially 
uniform density. It is straightforward to diagonalize the 
Hamiltonian and obtain the mean-field equations [9] 



Ck 



e 2 n 2 D 



C 



VA(k-k') 



1 - 



E' 



A k = - 



(3) 
(4) 



where £k — e k~ C = e/2ird is the capacitance per unit 
area, and #k,k' = $k ~ $k'- The form factor for the two 
stackings are 



/(0k,k' 



(1 + cos6» k ,k') 

COS #k,k' (1 + COS #k,k' 



Hexagonal , , 
) Bernal ' [ ' 



We point out that the self-energy in Eq.Q takes into 
account both intrinsic and extrinsic contributions that 
cancel the cos£?k,k'-dependent terms in the form factor 
and make the results independent of the ultra-violet cut- 
off mm. Therefore the intra-layer self-energy in Eq. ([3]) 
is the same as that for a conventional system [l9l . [20j | . The 
chemical potential /i is determined by the carrier density 
that takes into account the four-fold spin and valley de- 
generacy 



n 2 D 



1 



£k 



(6) 



It is straightforward to derive similar equations for a con- 
ventional electron- hole system [§] . They are obtained by 
changing the single-particle dispersion to a quadratic and 
replacing the form factor /(#k,k') by a constant / = 2. 
We solve Eqs. ©, (0| and © iteratively to obtain self- 
consistent results. 
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FIG. 1: (Color online) Excitonic gap Ak and the quasiparticle 
energy _E" k in graphene bilayer for Bernal (green dashed) and 
hexagonal (red solid) stacking with kpd — 1. The quasiparti- 
cle spectrum becomes linear with a renormalized velocity 
va > va for large k 2> &f- The dotted blue curves show 
corresponding results for an electron-hole system at r 3 = 2.7 
and kpd — 1 when plotted using relevant (atomic) unit for 
energy 0. 

Results: Figure [1] shows the excitonic gap Ak and the 
quasiparticle energy E k for the Bernal (green dashed) 
and the hexagonal (red solid) stacking. The excitonic gap 
Ak is maximum at the Fermi momentum kp where the 
quasiparticle energy E k is minimum. Since the electron- 
hole Coulomb interaction is always attractive, the exci- 
tonic condensate order parameter is nonzero down to the 
bottom of the Fermi sea, Ak=o 7^ 0. Our results predict 
that the hexagonal-stacked system will have a larger exci- 
tonic gap than the Bernal-stacked system. The quasipar- 
ticle energy E k becomes linear at large );> fcp, since the 
constituent particles of the exciton have a linear disper- 
sion. The speed of these quasiparticles is increased due 
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to intralayer exchange self-energy [20|, [2l|, |22| although 
the increase is modest, ~ 10%. Corresponding results 
for a conventional electron-hole system (blue dotted) are 
also shown in Fig. [TJ 
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FIG. 2: (Color online) Dependence of the graphene bilayer ex- 
citonic gap A m = Max(Ak) on interlayer distance d for Bernal 
(green dashed) and hexagonal (red solid) stacking. This gap 
can be tuned by changing nir> for a given sample. Corre- 
sponding result for a conventional system at r a = 2.7 is shown 
in dotted blue. 

Figure [5] shows the dependence of the maximum exci- 
tonic gap A m on kpd. We find that A m is weakly depen- 
dent on the stacking and decays rapidly when kpd^> 1. 
This result implies that the excitonic condensation is a 
robust phenomenon that will not require precise align- 
ment of the two graphene sheets when they are being em- 
bedded in a dielectric. With typical graphene carrier den- 
sities U2D ~ 10 12 /cm 2 and d ~ 100A or kpd ~ 1, the 
excitonic gap is appreciable, A m ~ 30 meV. 

A direct probe of the excitonic gap is the quasiparti- 
cle density of states. For graphene with no interactions, 
the density of states is linear, Dq(E) = 2E/ttH 2 Vq. In 
the excitonic condensate phase, for intermediate energies 
A m < E < i?k=o there are two rings in the phase-space 
consistent with that energy: one with k < kp and the 
other with k > kp. Therefore the quasiparticle density 
of states is given by D(E) = D < (E) + D > (E) where £>< 
(£>>) denotes the density of states from respective rings. 
Figure [3] shows D < (E) (green dashed) and D > (E) (red 
solid); they are both zero for E < A m and diverge at 
A m as is expected. Note that D < (E) = for E > i?k=o, 
since there are no states for k < kp with energies higher 
than -Ek=o- The asymmetry in D < and D > for E 3> A m 
is due to the linear dispersion of carriers and the nonzero 
electron-hole pairing that extends to the bottom of the 
Fermi sea, A k=0 ^ Q . The inset shows corresponding 
results a conventional system, where the density of states 
without interactions is constant, D (E) = m/ith 2 . 

Superfluidity of a uniform Bosc-Einstein condensate is 
characterized by a non-zero phase stiffness p s that quan- 
tifies the energy of a condensate with a linearly wind- 
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FIG. 3: (Color online) Quasiparticle density of states contri- 
butions D < (E) (green dashed) associated with states with 
k < kF, and D > (E) (red solid) associated with states k > kp. 
These results are for hexagonal stacked graphene bilayers with 
kpd = 1. Both diverge at E — A m , as expected. The total 
density of states D = D < + D > can be probed by differential 
conductance for tunneling from a metal into the condensate. 
The inset shows corresponding results for an electron-hole sys- 
tem at r s =2.7 and kpd = 1. All results are expressed in their 
respective units. 



ing phase, E(Q) = p s Q 2 A/2 where A is the area of 
the sample and the phase of the condensate varies as 
<&(x) — Qx. For graphene, since Ep is the sole energy 
scale (at zero temperature), it follows from dimensional 
analysis that phase stiffness must scale linearly with the 
Fermi energy, p s = g(kpd)Ep where g(x) is a dimension- 
less function that satisfies g ~ 0(1) [24| when < x < 1 
and g — * for x ^> 1. Hence, the phase stiffness is 
given by p s — g(kpd)hvc\/^n2D- The condensate en- 
ergy E(Q) can also be expressed, in the particle-picture, 
as the kinetic energy of excitons that have condensed 
in a state with center-of-mass momentum HQ. Thus, 
E(Q) = Nh 2 Q 2 /2m ex where m ex (N) is the mass (num- 
ber) of condensed excitons 



25j . Equating the two ex- 



pressions for energy implies m ex = n2Dh 2 / p s 
Thus we predict that the phase stiffness p s and the exci- 
ton mass will both vary as the square root of the carrier 
density. We emphasize that these results are unique to 
graphene and, as we will show in the next section, are 
equivalent to the absence of excitonic Wigner crystalliza- 
tion in graphene bilayers 2, 26|. 

Discussion: In this paper, we have investigated the prop- 
erties of excitonic condensates in graphene bilayers. Our 
calculations predict that excitonic condensation will oc- 
cur at all carrier densities as long as kpd ~ 1, and that 
the strength of the condensate, as measured by the exci- 
tonic gap A m is relatively insensitive to the stacking. 

The mean-field results presented in this paper are ob- 
tained at zero temperature T=0. (Finite temperature 
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analysis gives a critical temperature Tmf/Ef ~ 0.2 or 
Tmf ~ 20 meV. This is an artifact of the mean- field 
approximation.) In two dimensions, the critical temper- 
ature T c for Bosc-Einstein condensation is zero, but the 
superfluid properties survive for T < Tkt where Tkt is 
the Kosterlitz-Thouless transition temperature. There- 
fore, our results will be valid at nonzero temperature 
T -C Tkt [27|. A weak disorder will suppress the ex- 



was completed, we became aware of a recent related 
work 27 1 . 



citonic condensate order parameter and reduce the exci- 
tonic gap, an effect equivalent to increasing the value of 
kpd. Therefore, we have ignored the effects of a weak 
disorder potential. 

In our analysis, we have only considered excitonic con- 
densation with uniform density. In conventional (quan- 
tum Hall electron-hole) bilayers, varying d and r s (u) 
leads to excitonic condensates with lattice structure [lCj, 
I2H l28j |. The origin of the lattice structure is Wigner 
crystallization of carriers in an isolated layer at large r s 
(small v). Graphene does not undergo Wig ner crystal- 
lization as its carrier density is changed [29(. Therefore, 
we expect that the excitonic condensate in graphene bi- 
layers remains uniform. Now we show that this result is 
equivalent to our predictions for density dependence of 
p a and m ex . The quantum kinetic energy of an exciton, 
associated with localizing it within a distance 1/kp, is 
K = H 2 kp/2m ex . The potential energy due to the dipo- 
lar repulsion between them is P = e 2 d 2 k F /e. Hence their 
ratio is given by P/K = e 2 d 2 kFm ex /eh 2 . Wigner crys- 
tallization occurs when the ratio P/K ^S> 1. This ratio 
will solely be a function of kpd - no matter what the 
value of d is - if and only if m ex oc kp = yfrn^B. There- 
fore, results in the last section show that the excitonic 
condensate in graphene will not undergo Wigner crystal- 
lization in spite of the dipolar repulsion between excitons 
with a quadratic dispersion. This result, too, is unique to 
graphene and is markedly different from the behavior of 
dipolar excitonic condensates in conventional bilayers. It 
is interesting that the mass of these effective bosons has 
the same density dependence and order of magnitude as 
the cyclotron mass of fermionic carriers in graphene 

The onset of excitonic condensation can be detected 
by a divergent interlayer drag [30I ]. A uniform in-plane 
magnetic field Bu between the two graphene sheets is 
expected to induce a (counterflow) supercurrent Jd in 
such a condensate [3lj |. Jd — 2p s e 2 dBu/h 2 . The phase 
stiffness p s and its density dependence can be directly 
obtained from experimental measurements of the coun- 
terflow supercurrent. The verification (or falsification) 
of our predictions, including the density dependence of 
p s and m ex , will deepen our understanding of properties 
and condensation of excitons with massless fermions as 
constituent particles. 
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